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Abstract 

By using the method developed in the paper [Georg. Inter. J. Sci. Tech., Volume 3, 

Issue 1 (2011), 107-129], it is obtained a representation in an explicit form of the weak 
solution of a linear partial differential equation of the higher order in two variables with 
initial condition whose coefficients are real-valued simple step functions 
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1. Introduction 

In Q has been obtained a representation in an explicit form of the solution of the linear 
partial differential equation of the higher order in two variables with initial condition whose 
coefficients were real-valued coefficients. The aim of the present manuscript is resolve 
an analogous problem for a linear partial differential equation of the higher order in two 
variables with initial condition whose coefficients are real-valued simple step functions. 
The paper is organized as follows. 

In Section 2, we consider some auxiliary results obtained in the paper [5l|- In Section 3, 
it is obtained a representation in an explicit form of the weak solution of the partial differ¬ 
ential equation of the higher order in two variables with initial condition whose coefficients 
are real-valued simple step functions. 
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2. Some auxiliary results 

Definition 2.1. Fourier differential operator (^) in is defined as follows : 


0 0 0 0 0 0 0 

0 0 ^ 0 0 0 0 

0 -^ 0 0 0 0 0 

0 0 0 0 f 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 ^ 

00 0 0 0 0 

For n G N, let FD”[—/,/[ be a veetor spaee of all n-times differentiable funetions on 
[—/,/[ sueh that for arbitrary 0<^<n — 1, a series obtained by a differentiation term by 
term of the Fourier series of pointwise eonverges to for all x G [—/,/[. 

Lemma 2.2. Let f G FD^^'>[—l,l[. Let Gm be an embedding of the in to R°° 

which sends a function to a sequence of real numbers consisting from its Fourier coeffi¬ 
cients. i.e., if 



then Gf{f) = {^,c\,d\,C 2 ,d 2 ,...). Then, for f G [— the following equality 

holds. 

Proof. Assume that for / G FD^^'> [—l-,l[, we have the following representation 



By the definition of the elass FD^^'i [—/,/[, we have 
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By the seheme used in the proof of Lemma 2.1, we ean get the validity of the following 
assertion. 

Lemma 2.3. Let Gm be an embedding of the FD^ [~fl[ to R°° which sends a function to 

a sequence of real numbers consisting from its Fourier coefficients. 

Then, for f G FDG) [—/,/[ and Ak G R{0 < k < n), the following equality 

(G?‘o(i;A.((?)|-)‘)oG,)(/) = £»,^(/) (2.3) 

^k=(3 OX / / ox 


holds. 
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Example 2.4. ^ If A is the real matrix 


then 


e^^=e^^ 


la CO 
\-(o a 


(2.4) 

cos (cot) 

— sin(cot) 

sin(cot) \ 
cos (cot) J ' 

(2.5) 


Lemma 2.5. For m>\, let us consider a linear autonomous nonhomogeneous ordinary 
differential equations of the first order 


X {{ak)km) + {fk)keN ( 2 . 6 ) 

with initial condition 


{akiO))k(^fil = {Ck)keN, 

where 

(0 (Ci)it6N £ 

(//) / = {fk)ke'H ihe sequence of continuous functions of a parameter t on R. 
For each k>\,we put 


(2.7), 


m 


kn. 


^_ r\ ^ 


n—0 
m—l 


kn. 


^k=Yj (“l)”'^2«+l(y 


2n+\ 


n—0 


Then the solution of (2.6)-(2.7) is given by 


{ak{t))km = e 






( 2 . 8 ) 

(2.9) 

x/(x)^/T, (2.10) 


where e 




X {Ck)k€N Z e 


denotes an exponent of the matrix t {LZoAnimi)) and it 


exactly coincides with an infinite-dimensional (1,2,2,...) -cellular matrix D{t) with cells 
{Dk{t))ke'N for which Do{t) = {e*^°) and 


y—sin(c0i;f) cos(co^:f) J 


( 2 . 11 ) 


where for k>\, (5k and (Ok are defined by (2.8)-(2.9), respectively. 


Proof We know that if we have a linear autonomous inhomogeneous ordinary differential 
equations of the first order 


with initial condition 


-^^{{ak)k€n) = E X {{ak)keN) + ifk)keN 


{(lkiO))k(^fil = (Q)jt6N, 


( 2 . 12 ) 

(2.13) 


4 


where 

(i) {Ck)kefi G R”; 

(//) {fk))keN is the sequence of continuous functions of parameter t on 7?; 

{in) E is an infinite dimensional (1,2,2,... )-cellular matrix with cells {Ek)km- 
Then the solution of (2.6)-(2.7) is given by (cf. 01, §6, Section 1) 

{ak{t))ken = X {Ck)ken+[ x f{x)dx, (2-14) 

Jo 

where and denote exponents of matrices tE and (x — t)E, respectively. 

Note that t is an infinite-dimensional (1,2,2,... )-cellular matrix 

with cells {tEk)keN such that tEo = (tAo) and 


tEk={^^^ 1 (2.15) 

\-t(3ik tOk J ^ ’ 

for k > 1. Under notations (2.8)-(2.9), by using Example 2.4 we get that for t £ R, 
exactly coincides with an infinite-dimensional (1,2,2,...) -cellular matrix D{t) with cells 
{Dk{t))ken for which Do(f) = and 

D,(,) ). (2.16). 

^ ^ y—sm(cOi:t) cos{(i)kt) J ^ ^ 

Note that, for 0 < X < t, the matrix exactly coincides with an infinite-dimensional 

(1,2,2,...) -cellular matrix D{x — t). 

□ 


The following proposition is a simple consequence of Lemma 2.5. 

Corollary 2.6. Eor m > 1, let us consider a linear partial differential equation 

d 3” 

^'T(t,x) = £a„—'T( t,x) ((f,x) G [0,+oo[x[-/,/[) (2.17) 

at ax 

with initial condition 

'T(0,x) = y + ^^Qcos -h^sin e ED''^\-1,1[. (2.18) 

^(y ,ci,<ii,C 2 ,r/ 2 ,...) is such a sequence of real numbers that a series 'E(t,x) defined by 

+ (cicos((0)if)-h 

^ it=i ^ 

Jctzjc Jctzx \ 

(7i:Sin((0^t))cos(—^) -h {dkcos{o)kt) — Cksin(cOkt)) sin(—^)J (2.19) 

belongs to the class ED^^’^'>\—fl[ as a series of a variable x for all f > 0, and is differ¬ 
entiable term by term as a series of a variable t for all x £ [—l,l[, then 'T is a solution of 
(2.17)-(2.18). 
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3. Solution of a linear partial differential equation of the higher 
order in two variables with initial condition when coefficients 
are real-valued simple step functions 

Let 0 = <ti = T and —/ = xq < • • • < xy = /. Suppose that 


;=Oi=o 


For m > 1, let us consider a partial differential equation 


d 2m 

= Y^An{t,x)—^>{t,x) {{t,x) G [0,r[x[-Z,Z[) (3.1) 

at ax 


with initial condition 


'F(0,.r) = y+ ^c;tcos(^^^ +4sin(^^^ eFD^^\-l,l[. 
k—\ 


(3.2) 


Definition 3.1. We say that ^{t,x) is a weak solution of (3.1)-(3.2) if the following condi¬ 
tions hold: 

(i) *F(t,x) satisfies (3.1) for each (t,x) G [0,r[x[—Z,Z[ for which t / t,(0 < i < I) or 

X / Xj{0 < 7 < /); 

(ii) *F(t,x) satisfies (3.2); 

(iii) for each fixed x G [—Z,Z[, fhe funcfion 'F(t,x) is continuous with respect to t G[o,r[, 
and for each t G [0, r[ the function *F(t,x) is continuous with respect to x on [—Z,Z[ except 
points {xy : 0 < 7 < 7 — 1}. 

First, let fix j and consider a partial differential equation 

d I'n ^ d” 

-'F(oj)(z,x) = lAf^^(z,x)^'F(oj)(z,x) ((t,x) G [0,+oo[x[-Z,Z[) (O.MPDE) 
with initial condition 


'^{o.j){to,x) = j + ^ckcos(^^'^+dksmy ^ j 


k=i 


/kTlX\ _ 


JOd) 


. f;4»J'cos(f^) +4»''sin(fH) £(oj)(;c) 

k—l 


2 I ) 

By Corollary 2.6, under some restrictions on {^,c\,d\,C 2 ,d 2 ,- ■ ■), a series *F(oj)(z,x) 
defined by 

tA^j) (0./) 


^(0,j){t,x) = 


k=l 


df'^^ sin(co^°’'^V))cos(^^) -|- {df’^^co&{(i)f’^\) — sin(co^°’''\)) sin(^^)^ (3.3) 
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is a solution of {0.j){PDE)-{0,j){IC), 

Now let consider a partial differential equation 


iit,x) € [0,+oo[x[-Z,Z[) (I.MPDE) 


with initial condition 


We wiil try to present the solution of the {\,j){PDE) by the following form 




(ij) 

e“o c. 






k=\ 


sin(co^^’-^^f))cos(^^) + (rZj[^’‘^^ cos(co[^’'^\) — sin(co^^’''^t))sin(^^)^ (3.4) 

In order to get validity of the condition (l,y)(/C), we consider the following infinite 
system of equations: 




(3.5) 


e®* cos(co[^’'^\i) + sin((o['’'^\i)) = 

cos((of’^'^Zi) +4°’^'^ sin((of ^')zi))()t G N), 


(3.6) 


(4^’'^'^ cos(co4''4) — sin(co4'^4)) =e'^* (4°'^^ cos(co^°’'^4) — 4°’'^^ sin(co^°’''4))(^ G N). 

(3.7) 

We have 

/i . . (O.i) . (1.;K /n .•^ 

(3.8) 


JiJ)_„n(Ar->-A<‘-)) (OJ) 
Cq — e ^0 • 


For ^ G N we can rewrite equations (3.6)-(3.7) as follows: 


cos(co4-^4) +4^’'^'^sin((o4-^4)) = (4°'^'^cos((o[°’-^'4) +4*^’'^^sin(co4-^'4)), 


(3.9) 






—4^ sin ((o4''4) + 4^ cos (co^-^^ Zi) = 

A = (4°’'^'^cos((o[°’‘^^Zi) +4°’'^'^ sin(co[°’-^'4)) (3-11) 

and 

B = (4°’^'^ cos((of ^'4) - 4°’^^ sin((oPzi)): 


Setting 


)f’-'^Zi)-4°’^4in((o[°’^^Zi)). 

(3.10) 


(3.12) 
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for ^ G N we obtain 



eos(co[^’-^'^Zi) + sin(co[^’-^'^fi) = A 

(3.13) 

and 

—c^^ sin(co[^ 6) + eos (co[^ ) = B • 

(3.14) 

It is obvious that the system of equations (3.13)-(3.14) has the unique solution whieh ean 
be done as follows: 

= Aeos(co[^’'^'^Zi) —Bsin(co[^’-^^C) (3.15) 

and 

= ]Bcos(co[^’'^'^Zi) + Asin(co[^’-^'^fi) 

(3.16) 

for Z: G N. 

By Corollary 2.6, under some restrietions on {%^,c^i’'’\d[^’^\c 2 ’'’\d 2 ’^\- ■ 

,.), the se- 


ries *P(i defined by (3.4) is the solution of {\.j){PDE)-{\,j){IC), 

It is obvious that under niee restrietions on eoeftieients partieipated in (3.1) and (3.2), 
we ean eontinue our proeedure step by step. Correspondingly we ean eonstruet a sequenee 
(*P (5 y))o<i</-i,i<j< 7 -i sueh that satisfies a linear partial differential equation 


3 ^ {k,j), . 3 ' 


^'F(,,,.)(t,x) = IAr’^^(t,x)^'F(,,,.)(t,x) ((t,v) G [0,+oo[x[-Z,Z[) {s,i){PDE) 


n=0 


with initial eondition 


'3x' 






k=i 


V Z J 


Now it is obvious to observe that we have proved the validity of the following assertion. 

Theorem 3.2. If for coefficients (^, C 2 ,d 2 '^\.. .){\ <i <f\ < j <J) func¬ 
tions {fx) satisfy conditions of Corollary 2.6, then a function ,x) : [0, r[x[-z,zK 
R defined by 

i-ij-i 

II '^{ij)ixf) '>^lL[ti,ti+i[x[xj,Xj+i[ifiX) (3-17) 

;=0 ;=0 

is a weak solution of (3.1 )-(3.2). 


Example 3.3. Let consider a linear partial differential equation of the 22 order in two 
variables 

3 3^ 3^^ 

^*E(z,v) =A{t,x) X ^*E(t,x) +B{t,x) X ^-^*E(t,v) ((z,x) G [0,27t[x [0,7t[) (3.18) 

with initial condition 

*E(0,x) = +5sin(v), (3-19) 
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where 

A[t,x) = X[o,7i[x[0,7t[(^)''') A 1.55 X X[Ti^2n[x[0,n[{f yX) 

and 


B{t ,x) — 2 X X[0,7l[x[0,7t[(l)''') 2 ^ X[7t,27t[x[0,Jt[(l )''-)• 

The programm in MathLab for a solution o/(3.18) — (3.19), has the following form: 

A1 = [0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 

A2 = [0,1.55,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2]; 

Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 

D1 = [5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

AlO = 0;A20 = 0;C10 = 0.015; 
for ^=1:20 
S\{k) =A\0-,S2{k) =A20; 
for « = 1 : 10 

51 (^) = 51 (^) + (—1)U) *Al(2*?i) *k^2*n); 

S2{k) = S2{k) + (— 1)U) *A2{2*n) *k^2*n)', 

end 

end 

for ^=1:20 
O\{k)=0- 
02{k) = 0; 

end 

for ^=1:20 
for « = 1 : 10 

Ol{k) = Ol{k) + {-lY*Al{2*n + l)*k^2*n + l)\ 

02lk) = 02{k) + {-lY*A2{2*n + Y*k^2*n + Y\ 

end 

end 

[Tl,Xl\ =meshgrid(0 : {pi/lO) : pi,0 : {pi/\0) : pi); 

Zl = 0.5*C10*exp{Tl.*Al0); 
for ^=1:20 

Z\=Z\+C\{k)*exp{T\*S\{k)).*cos{X\.*k).*cos{T\*0\{k))+D\{\)*exp{T\* 
51 (^)). >i:cos(Xl. *k).* sin(ri * 0\{k))+ 

D\{k) >i:exp(ri >i:51(^)). * sin(Xl. *k).* cos(ri * 01(^)) — Cl(^) * exp(ri * 51(^)). * 
sin(Xl. *k).* sin(ri * 01(^)); 

end 

C20 = exp (p/ * (A 10 - A20)) * CIO; 
for ^=1:20 

A(^) = exp((51(^) — S2{k)) * pi) * {C\{k) *cos(d(^) * pi) +Dl(fc) >i:sin(01(^) * pi)); 
B{k) = exp((51(^) — S2{k)) * pi) * (Dl(^) *cos(d(^) * pi) — Cl(^) *sin(d(^) * pi)); 

end 

for ^=1:20 

C2(^) =A{k) *cos{02{k) * pi) — B{k) * sm{02{k) * pi); 
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D2{k) = B{k) *cos{02{k) * pi) +A{k) * sm{02{k) * pi)-, 

end 

[T2,X2] =meshgrid(p/: (p//10) : {2*pi),0 : (p//10) : pi)-, 

Z2 = 0.5 *C20 * exp{{T2) *A20); 
for ^=1:20 

Z2 = Z2+C2{k)*exp{T2*S2{k)).*cos{X2.*k).*cos{T2*02{k)) +D2(1) =t:exp(r2=t= 
SI (A:)). =t=cos(S'2. *k).* sin(r2 =t= 02{k))+ 

D2{k) *e-xp{T2*S2{k)). * sin(S'2. *k).* cos(r2 * 02{k)) — C2{k) * exp(r2 * S2{k)). * 
sin(X2. *k).* sin(r2 * 02{k))-, 

end 

surf(ri,Xl,Zl) 

hold on 

surf(r2,Z2,Z2) 

hold off 


Picture 1 


6 

5 

4 

N 3 

2 

1 

0 

4 
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1. Graphic of the solution of the LPDE-(3.18) with IC-(3.19)). 


Example 3.4. Let consider a linear partial differential equation of the 22 order in two 
variables 


dt 


*P(r,x) = A(f,x)*P(i,x)+B(i,x) X 


3x2 


'P(r,x) +C(i,x) X 


3x2 


'P(r,x)+ 


10 











Picture 2 
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4 
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-6 

4 
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2. Graphic of the solution of the LPDE-(3.20) with IC-(3.21)). 




with initial condition 


D{t,x) X {{t,x) G [0,27i[x[0,7t[) 


*^(0,^:) = +5sin(x), 


where 


and 


A{t,x) = lX[0,7c[x[0,7t[(l)-^) 3“ 0 X 2 ji:[x [0,7c[(E-^)) 
B{t,x) = X[0,7t[x [0,71[(E''') + 0 X X[jl,27l[x[0,7t[(^)''-)) 
C(t,x) = 0 X X[0,7t[x [0,7t[(E''') + 1 X X[7C,27t[x [0,Jt[(l)''') 

D(t,x) = 2 X X[o,jt[x [o,7t[(E''') + 2 X X[ 7 t, 27 t[x [o,7i[(^)''') 


(3.20) 

(3.21) 


The graphical solution of(3.20)-(3.21) can be obtained by MathLab programm used in 
Example 3.3 for the following data: 

A1 = [0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 

A2 = [0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 

Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 
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D1 = [5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 
AlO = 1;A20 = 0;C10 = 0.015; 


Example 3.5. Let consider a linear partial differential equation of the 21 order in two 
variables with constant coefficients 

=A{t,xmt,x) + ^^>{t,x) + l^^>{t,x) {{t,x) G [0,27i[x[0,7i[) (3.22) 

with initial condition 


Picture 3 


150 

100 

50 

M 0 

-50 


-100 

-150 
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8 



3. Graphic of the solution of the LPDE-(3.22) with IC-(3.23)). 

*P(0,x) = -|-5sin(.r), (3.23) 

where 

A(t,x) = lX[0,7l[x[0,Jt[(l)''') A 0 ^ X[jt,27C[x 

Since 

1 = 1 X X[0,7C[x [0,7C[(E-^) + 1 X X[7t,27t[x [0,7C[(l a) 

and 

2 = 2 X X[0,7t[x [0,7c[(l)''') + 2 X X[7C,27C[x [0,7t[(l)-^)) 

the graphical solution of (3.22)-(3.23) can be obtained by MathLab programm used in 
Example 3.3 for the following data: 
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A1 = [0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 
A2 = [0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 
Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

D1 = [5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 
A10= 1;A20 = 0;C10 = 0.15; 


Example 3.6. Let consider a linear partial differential equation of the 21 order in two 
variables 


a 

dt 


'E(r,x) = A{t,x)^{t,x) +B{t,x) X 


a2 




a2i 




'E(cx) ((Cx)g[0,27i[x[0,7i[) 
(3.24) 


Picture 4 


N 


4. Graphic of the solution of the LPDE-(3.24) with lC-(3.25)). 


with initial condition 


*P(0,.r) = ^ + 100sin(:t), 


(3.25) 


where 

and 


A(t,x) — lX[0,7[[x [0,71:[(E-^) A 0%[;j_27t[x [0,7t[(C-^) 
B{t,x) = X[0,7C[x[0,7t[(C-^) ~ X[7C,27 C[x [0,7C[(C''') ■ 
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The graphical solution of(3.24)-(3.25) can be obtained by MathLab programm used in 
Example 3.3 for the following data: 

A1 = [0,1,100,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 

A2 = [0,-1,100,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 

Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

D\ = [100,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

A10= 1;A20 = 0;C10 = 0.15; 

We see that we have no graphic on the region [ti, 271 [x [0,71 [ which hints us that coeffi¬ 
cients of the LPDE (3.24)-(3.25) on that region do not satisfy conditions of Theorem 3.2. 

Example 3.7. Let consider a linear partial differential equation of the 21 order in two 
variables 

^W{t,x)=A{xf)W{t,x)+50^W{t,x) + ^W{t,x)+2^W{t,x){{t,x)e[0MA^A) 

(3.26) 

with initial condition 


Picture 5 



5. Graphic of the solution of the LPDE-(3.26) with lC-(3.27)). 

*P(0,x) = +5sin(.r), (3.27). 

where 

A(t.,x) = X[ 0 ,n:[x[0,7t[(E-^) +0 X X[7i:.27t[x[0,7t[(C''')- 
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The graphical solution of(3.26)-(3.27) can be obtained by MathLab programm used in 
Example 3.3 for the following data: 

A1 = [0,0,0,0,50,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 

A2 = [0,0,0,0,50,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0]; 

Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

D\ = [7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

A10= 1;A20 = 0;C10 = 0.15; 


Picture 6 



X 0 0 


6. Graphic of the solution of the LPDE-(3.28) with lC-(3.29)). 

Example 3.8. Let consider a linear partial differential equation of the 22 order in two 
variables 

d d d^^ 

=A{x,t)'¥{t,x)-^'¥{t,x)+2^'¥{t,x) {it,x) G [0,27t[x[0,7i[) (3.28) 

with initial condition 

*P(0,;r) = +5sin(x), (3.29). 

where 

Af,x) = X[0,7t[x [0,71[(E''') +0 X X[jI,27l[x[0,7t[(C'^)- 

The graphical solution of(3.28)-(3.29) can be obtained by MathLab programm used in 
Example 3.3 for the following data: 
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A1 = [-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 
A2 = [-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 
Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,01; 

D1 = [5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

AlO = 1;A20 = 0;C10 = 0.015; 


Remark 3.9. Notice that for each natural number M > 1, one can easily modify the MathLab 
program described in Example 3.3 for obtaining the graphical solution of the linear partial 
differential equation (3.1)-(3.2) whose coefficients (A„(t,x))o<H< 2 M are real-valued simple 
step functions on [0, r[x [—/, /[ and / is a trigonometric polynomial on 
Remark 3.10. Since each constant c same times is a step function we can use MathLab 
program described in Example 3.3 for a solution of the linear partial differential equation 
(2.17)-(2.18) with constant coefficients. 

Eor example, if we consider LPDE 

= + ((t,x) E [0,27i[x [0,7t[) 

^ ^ ^ ^ ^ (3.30) 


Picture 7 


M 


7. Graphic of the solution of the LPDE-(3.30) with 1C-(3.1)). 
with initial condition 

*P(0,x) = + 150sin(x), 


(3.31) 
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then for a solution (3.30)-(3.31),in MathLab programm described in Example 3.3 we must 
enter the following data: 

A1 = [-1,0,1,0,-1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 

A2 = [-1,0,1,0,-1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2]; 

Cl = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

D\ = [150,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

AlO = 1;A20 = 0;C10 = 0.015; 

Remark 3.11. The approach used for a solution of (3.1)-(3.2) can be used in such a case 
when coefficients (A„(t,x))o<«< 2 M are rather smooth continuous functions on [0, T [x[-Z,Z[. 
If we will approximate (A„(t,x))o<H< 2 M by real-valued simple step functions, then it is nat¬ 
ural to wait that under some ’’nice restrictions” on (A„(t,x))o<«< 2 M the solutions obtained 
by Theorem 3.2, will give us a ’’good approximation” of the solution of the required linear 
partial differential equation of the higher order in two variables with corresponding initial 
conditions. 
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